De Sitter Universe in Nonlocal Gravity 



E. Elizalde^*, E.O. Pozdeeva^t, and S.Yu. Vernov^'^^ 

-'^Instituto de Ciencias del Espacio (ICE/CSIC) and 
Institut d'Estudis Espacials de Catalunya (lEEC), 
Campus UAB, Facultat de Ciencies, Torre C5-Parell-2a planta, 
E-08193, Bellaterra (Barcelona), Spain 
^Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, 
Leninskie Gory 1, 119991, Moscow, Russia 

- ■ Abstract 

A nonlocal gravity model, which does not assume the existence of a new dimensional 
parameter in the action and includes a function /(D^^i?), with □ the d'Alembertian 
operator, is considered. The model is proven to have de Sitter solutions only if the 
function / satisfies a certain second-order linear differential equation. The de Sitter 
solutions corresponding to the simplest case, an exponential function /, are explored, 
without any restrictions on the parameters. If the value of the Hubble parameter 
is positive, the de Sitter solution is stable at late times, both for negative and for 
positive values of the cosmological constant. Also, the stability of the solutions with 
zero cosmological constant is discussed and sufficient conditions for it are established 
in this case. New de Sitter solutions are obtained which correspond to the model 
with dark matter and stability is proven in this situation for nonzero values of the 
cosmological constant. 

1 Introduction 

Modern cosmological observations, such as those coming from Supernovae la (SNe la) [1], 
from the cosmic microwave background (CMB) radiation [21 13], large scale structure (LSS) 
baryon acoustic oscillations (BAO) [5], and weak lensing [6], allow to obtain joint constraints 
on cosmological parameters (see, for example, [7]) and indicate that the current expansion of 
the Universe is accelerating. There are a few types of models able to reproduce this late-time 
cosmic acceleration. The simplest one is general relativity with a cosmological constant (for 
a review, see, e.g., [8]). Some others involve modifications of gravity, as for instance F{R) 
gravity, with F{R) an (in principle) arbitrary function of the scalar curvature R (for recent 
reviews see, e.g., [9| flO]). 
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Modified gravity cosmological models have been proposed in the hope of finding solutions 
to the important open problems of the standard cosmological model. There are lots of 
ways to deviate from Einstein's gravity. Different modifications of this theory have been 
considered in detail in the reviews [HI [TT]. As a promising modification of gravity, the 
nonlocal gravity theory obtained by taking into account quantum effects has been proposed 
in [12]. Also, as is well known string/M-theory is usually considered as a possible theory 
for all fundamental interactions, including gravity; again, the appearance of nonlocality 
within string field theory is a good motivation for studying nonlocal cosmological models. 
Moreover, there was a proposal on the solution of the cosmological constant problem by a 
nonlocal modification of gravity [13]. The majority of nonlocal cosmological models explicitly 
include a function of the d'Alembert operator, □, and either define a nonlocal modified 
gravity [Tl[T5l[T6l[I7l[l8l|20l[l9l|2ll|22]oradda nonlocal scalar field, minimally coupled 
to gravity [23]. 

In this paper, we consider a modification that includes a function of the D"^ operator. 
Such modification does not assume the existence of a new dimensional parameter in the ac- 
tion. This nonlocal model has a local scalar-tensor formulation. A modification of nonlocal 
gravity with a term has been studied in order to realize a unified scenario, com- 

prising both early-time inflation and the late-time cosmic acceleration. It has been shown 
in [15] that a theory of this kind, being consistent with Solar System tests, may actually 
lead to the known Universe history sequence: inflation, radiation/matter dominance and a 
dark epoch. An explicit mechanism to screen the cosmological constant in nonlocal gravity 
was discussed in [TH [201 EI] • Different cosmological aspects of such nonlocal gravity models 
have been studied in[T6l[I71[T8l[20l[l9l[2ll[22], too. 

In this paper, we explore in detail de Sitter solutions in the nonlocal gravity model. We 
prove that the model can have de Sitter solutions only if the function f[n~^R) satisfies a 
certain second order linear differential equation, which is a nice result. The simplest and 
most studied model [151 US EZl UHl |20l HH, [21] admitting de Sitter solutions is characterized 
by a function f(n~^R) = f^e^^ where /o and f3 are real parameters. A few de Sitter 

solutions for this model have been found in [T5] and also analyzed in [20j. In both papers 
the authors put strict restrictions on arbitrary parameters (integration constants). 

Here, we will consider all possibilities for de Sitter solutions without any restriction. We 
will also obtain de Sitter solutions in the case when the matter included in the model is 
dark matter. Finally, we will analyze the stability of these de Sitter solutions and get the 
corresponding restrictions on the parameters of the model. In the case A = 0, the system of 
equation, describing this model, has been written in terms of Hubble-normalized variables 
and the stability of the fixed point of this system has been analyzed in [12]. This has been 
continued in [16] (see also [E]). We will consider this specific case in Subsection 14. 2[ 

The paper is organized as follows. In Section 2, we shortly review nonlocal gravity models 
that have a local scalar-tensor formulation. In Section 3 we obtain the necessary conditions 
on the function /(□~^i?), which allow us to get de Sitter solutions, and then look for general 
de Sitter solutions in the case of the exponential function f(n~^R). In Section 4 we discuss 
the stability of the de Sitter solutions for Friedmann-Lemaitre-Robertson-Walker (FLRW) 
and Bianchi I metrics. Section 5 is devoted to conclusions. 
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2 Nonlocal gravitational model 

Consider the following action for nonlocal gravity 



S = j rf^xv^ 1^ [R (1 + fia-'R)) - 2A] + £„| 



(1) 



where = Stt/Mpi^ the Planck mass being Mpi = G^^/^ = 1.2 x 10^^ GeV. We use 
the signature (—,+,+,+), g is the determinant of the metric tensor g^^, f a differentiable 
function, A the cosmological constant, and ^matter is the matter Lagrangian. Recall the 
covariant d'Alembertian for a scalar field, which reads 

D^^d,{V^g^^d,). (2) 
Introducing two scalar fields, t] and ^, we can rewrite action ([1]) in the following local form: 



S = j rf^a;v^|^[i?(l + /(77)-0+en^-2A] + /:^attcr| 



(3) 



By varying the action (j3]) over ^, we get □?] = R. Substituting rj = O ^R into action (|3]), 
one reobtains action ([T]). Varying action with respect to the metric tensor g^^, one gets 

Ig,, [R (1 + fir]) - - dp^d^v - 2A] - R^, (1 + /(r/) - + 

2 1 (4) 

+ 2 (^M^^-^ + ^M^^^^O - i9^.ua - V^d,) ifiri) - + /^'T^atter^. = , 

where is the covariant derivative and Tmatter^ii^ the energy-momentum tensor of matter, 
defined as 

rrj 2 -^matter 



- matter /lu 



(5) 



On the other hand, variation of action ([3]) with respect to rj yields 

□e + /'(r/)i? = 0, (6) 

where the prime denotes derivative with respect to rj. 
We take the spatially flat FLRW metric, 

ds'^ = -dt^ + a^it) [dxj + dxl + dxl) , (7) 

and consider the case where the scalar fields r] and ^ depend on time only. In the FLRW 
metric, system of Eqs. (j4]) is equivalent to the following equations: 

- 3i7^ (1 + /(r^) - + - ^Hj^ ifir]) - + A + /^Vm = , (8) 
(2H + (1 + fiv) -0 + + + 2i^|) ifiv) - - A + '^'^m = , (9) 



where H = a /a is the Hubble parameter, the dot denoting time derivative. For a perfect 
matter fluid, we have Tmattcroo = Pm and Tmatterii = Pradij- The equation of state (EoS) is 

p^= -3H{P^ + p^). (10) 

Adding up Eqs. ([8]) and IQ, we get 

2H (1 + fiv) - + + - Hj^ ifiv) - + ^'iP.n + pj = 0. (11) 

Furthermore, from Drj = R and ([6]) the equations of motion for the scalar fields rj and ^ 
follow 

i] + 3Hfi = - 6 f if + 2H^) , (12) 



^ + 3H^ = Q(^H + 2H''^f{r]), (13) 

where we have used R = 6H + 12H^. 

Let us consider the system of equations f|T0l) - f|T3|l . Together with ([8]), they are equiv- 
alent to the full system of Einstein's equations. Differentiating ([8]) with respect to t and 
substituting into Q), ffTOj) . f|T2|) . and f|T3|) . we get that ([8]) is an integral of motion for the 
system of equations ffTUl) - f[T51) . Therefore, to find a solution of the Einstein equation one can 
solve the system f[TU]) - f[T^ . which does not include the cosmological constant A. After that, 
substituting into ([8]) the initial values of the solution obtained, one gets the corresponding 
value of A. 

The system of equations considered does not include the function rj, but only /(r^), /'(^) 
and time derivatives of 1]. This property can be used to analyze the stability of the de Sitter 
solutions. 



3 De Sitter solutions 

3.1 Nonlocal models with de Sitter solutions 

We now assume that the Hubble parameter is a nonzero constant: H = Hq. In this case, 
Eq. ( fT2l) has the following general solution: 

7]{t) = - AHo{t - to) - r7oe-3^«(*-*°), (14) 

with integration constants to and rjo. All equations are homogeneous, so if a de Sitter solution 
exists at to = 0, then it exists at an arbitrary to- So, without loss of generality we can set 
to = 0. 

Subtracting Eq. ([8]) from Eq. IQ, we get a linear differential equation 

^ + 5Ho^ + QHlil + ^) - 2A + k\w,,, - l)pr^ = 0, (15) 
where ^(t) = f{ri{t)) - ^(t). 
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Solving (fT5|) and substituting ^(t) = fijjit)) — '^{t) into Eq. (fT3|) . we get a linear differ- 
ential equation for f{rj) 

V'f'iv) + + ^Hov - l2Hl) f'{r]) = ^ + SH^^ . (16) 

Therefore, the model, which is described by action (jS]), can have de Sitter solutions only if 
f{rj) satisfies Eq. (fT6|) . In other words Eq. ( [T6|) is a necessary condition that the model has 
de Sitter solutions. To prove the existence of de Sitter solutions for the given /(r^) one should 
also check Eqs. ([8]) and ([9]). Note that Eq. ( !T6|) has been obtained without any restrictions 
on solutions and the perfect matter fluid. 

To demonstrate how one can get /(//), which admits the existence of de Sitter solutions, 
in the explicit form, we restrict ourselves to the case ?7o = 0. In this case, Eq. ( IT6l) has the 
following form: 

16H'j"{v)-2AH'of'iv) = m, (17) 
where $(■?]) = ^{-AHot) = ^ + SHq^. We get the following solution 

fiv) = / |/ me-'^^'dC + leCsi/o'l e'^/'dC + C„ (18) 

where C3 and C4 are arbitrary constants. We can fix C4 without loss of generality. Indeed, 
it is easy to see that we can add a constant to f{ri) and the same constant to ^, without 
changing of Eqs. (jlj) and (E]). 

Following [20], we consider the matter with the EoS parameter w^a = -Pm/Pm to be a 
constant, not equal to —1. Thus, Eq. (|T0|) has the following general solution 

Pm = Poe-3(^+'"-)^°*, (19) 

where po is an arbitrary constant. 

Equation ( JTSll has the following general solution: 



*i(t) = Cie-^^^"* + C2e-2^°* - 1 + (20) 



At Po = 0, 

A 



At Po 7^ and Wm = 0, 

v]>2(t) = Cie-3^°* + C2e-2^"* -^ + ^2- ^e-3^"*t, (21) 



3Hn Hn 



,2 

At Po 7^ and = —1/3, 



Mt) = Cie-3^°* + ^26-2^0* - 1 + TT^I + (22) 
At Po 7^ 0, Will 7^ and Wm 7^ —1/3, 
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where Ci and C2 are arbitrary constants. 

Substituting the exphcit form of ^'(t), we get 

• For the model without matter (po = 0, = 

fiiv) = Y'"^' + Cse''^' + C„ (24) 

where C3 and C4 are arbitrary constants. Note that C2 is an arbitrary constant as 
welL 

• For the model with the dark matter {wm = 0, \E'(t) = \E'2(t)), 

f2{v) = m - '^e^'^^ (25) 

• For the model, including the matter with w^i = —1/3 (^(^) = ^3('^))) 

/3(^) = m + ^(^-lv] e^"- (26) 



• For the model, including the matter with another value of {^{i) = ^4(i)), 

fM = fM - ^^.f^.rr. e'^'"-^'^'''. (27) 

3(1 + 3Wra)H^ 

One can see that the key ingredient of all functions fiirj) is an exponent function. 

3.2 De Sitter solutions for exponential f{r]) 

In the following we consider de Sitter solutions for the model with 

fiv) = he^'^ , (28) 

where /o and /3 are constant. We choose this form of /(r^) because it is the simplest function, 
which belongs to the general set of functions fl2^ . fl25l) . and fl27|) . described above. Note that 
models involving this exponential form of f{ri) are being actively studied in the literature [151 

[isiiniiisiEniiisiEi]. 

In the case A = 0, the corresponding system of equation can been written in terms of 
Hubble-normalized variables. The stability of the fixed point for this system has been dis- 
cussed in [15] and the stability analysis for the same case has been further investigated in [16] 
(see also [I9]). We will consider this specific case in Subsection 14.21 In absence of matter, 
expanding universe solutions a oc have been found in [TJl |2T] . In [T7] the ensuing cosmol- 
ogy at the four basic epochs: radiation dominated, matter dominated, accelerating, and a 
general scaling has been studied for an interesting nonlocal model involving, in particular, 
an exponential form of f{ri). Screening of the cosmological constant in the nonlocal model 
described by the action ^ with an exponential f{ri) has been considered in [181 
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In [20], the model has been considered in further detail: the corresponding de Sitter 
solutions have been obtained, a condition to avoid a ghost has been considered, and a 
screening scenario for a cosmological constant has been discussed. Note that in [151 EO] 
the authors put restrictions on arbitrary parameters, to get de Sitter solution. It will be 
interesting to get de Sitter solutions without any restriction. 

Substituting the solution f lT^ of Eq. f|T2|) into Eq. f|T3|) and assuming that the integration 
constant r/o 7^ 0, we obtain that Eq. f lT3|) has the general solution 



(29) 



^(t) = 12^^^ [ i Ci + /"e(-'"'''''P[-^^o(*2-to)]-4//o{t2-to))//3+3jyot2^^2 j e-3^o*Mti 
with Ci and .^o arbitrary constants. If /3 = 2/3, then ^(t) can be found explicitly: 



The solutions obtained include four arbitrary parameters, namely rjo, ^o, Ci, and to- As we 
have mentioned above one can set to = 0, without loss of generality. 
At H{t) = Ho, Eq. dUD has the form 

+ ^fiv) (^f + v]- ^Hofiv)v - e + Hoi + + t^m)Pm = , (31) 



and, using f lT3|l . 

(f, + AHo)i + ^ + ri-Hofi- 12H^j + wJp^ = . (32) 

Substituting the explicit expression for t], we get 

3i/or/oe-'^"*e + ^fiv) (g^'e"*'^"* - 12r]oW + 2)6"=^^°* - 8(/3 - 2)) + /^^(l + w,^)p^ = , 
where 

From ([29D, it follows that 

^ = I2^^ll^e'^^°^ + y e[-''o'=''P(-^^o*2)-4Hot2]//3+3Hot2 ^^2^ _ (33) 

A straightforward calculation shows that Eq. (132|) has no solution for any value of the 
parameters such that fo ^ 0, rjo 0, and Hq ^ 0. 

Therefore, without loss of generality we can put r^o = 0. In this case, for /3 ^ 4/3, from 
(fT2l) and (fT3l) the following solution is obtained [2QJ: 

i = - ^e-'^^^^-'^yP + ^e-^^o(*-o) - eo, = - 4i/o(t - to), (34) 
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where Cq is an arbitrary constant, 

A = 3i7o(l + -^o), Po = ' " ^ ^ 3^' ^ ^ 

The case (3 = 2 corresponds to po = 0. Note that /3 = 2 corresponds to = —1/3. It 
means that the model with exponential f{ri) has no de Sitter solutions if we add this kind 
of matter. The type of function firf) that can have such solutions is given by (!26|) . 

For [3 = 4/3, we get 

m = - /o(co + 3Hoit - to))e-3^«(*-*°) - ^o, (36) 

A = 3H',{1 + eo), = 0, p^= - lff„Voe-^^°(*-*»). (37) 

This solution clearly corresponds to dark matter, because = 0. 



4 Stability of the de Sitter background 

4.1 The case of nonzero A 
4.1.1 The FLRW metric 

Let us now introduce new variables 

= fiv) = /oe''/^ iJ = V- (38) 
The functions 0(t) and iplt) are connected by the equation 

= ^<P^- (39) 

The system (|TTl) - f|T3l) can be expressed as a system of first-order differential equations, in 
terms of new variables. We rewrite Eqs. f lT^ and f[T^ as 



^P= - SH^p -6\^H + 2H'j , (40) 
i = ^, § = - ?,H§+^{h + 2H^^(P. (41) 

Using 

we get that Eq. ffTT|) is equivalent to 



2H (^1 + ^0 - = 4i7 - - -<pr + jH'<P -H- + ^m)Pm. (42) 
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Consider the de Sitter solution 



Pn 

For /3 7^ 4/3, we have 



and, for /3 = 4/3, 
As t tends to +oo. 



3(1 + 



(43) 
(44) 



3/o/3 4Hp(t_to)//3 _Co_ 3H„(t-to) _ f 

3/3-4 +3i/o ^' 



/o(co + 3Ho{t - to))e 



-3Ho(t-to) 



(45) 



for all -f^o > and /3 > 0. This system has a fixed point: 

= 0, e= -eo, ^= -4/^0, Pm = o. 

Note that we cannot fix Hq, using the relation A = 3Hq{1 + ^o), since is an arbitrary 
parameter. So, we have no isolated fixed point. 

Two different cases appear: A = and A 7^ 0. In this subsection, we consider the case 
A 7^ 0. For A = 0, the stability can be analyzed using a change of variables. This analysis 
will be presented in the next subsection. 

For A 7^ 0, one gets ^0 7^ in the neighborhood of the fixed point 



1 + 



13-Q 



1 + ^0^0 



(46) 



and we can divide Eq. (1421) in this expression to get the equation in the standard form: 

1 



H 



2 (1 + 



AH 



/3 



(47) 



We consider the time domain {ti, +00) such that l + — ^ 7^ 0. This expression is positive 
for A > and negative for A < 0. In the neighborhood of the fixed point, which corresponds 
to de Sitter solution, we have 



H{t) = Ho + 6h{t) + Oie^), 

0(t) = eMt) + 0{e^), 

ij{t) = - AHo + eMt) + 0{e^) 

Pm{t) = 6pmi{t) + 0{e^), 



(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
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where e is a small parameter. 

From Eqs. ([39]), (iO]), (jH]), and we obtain, to first order in e, the following 

system: 

4 

Pml = - ^-f^oPml, (54) 

01 = - 1^001, (55) 
= - 3//o^i + ^i/o'0i, (56) 



;i+eo 



12 

= - 3ifo^i - 12ifo/ii 



2 . 2 \ 



(57) 



(5J 



Note that the function is not included in this system. It can be defined using Eq. ([H]). It 
is plain that cannot tend to infinity, if all other first-order corrections are bounded. 
Let us now consider the system (!5^ - (!58l) . The functions 

Pmiit) = rfoe-^^°*/^ Mt) = rfie-^^«*/^ (59) 

where do and di are arbitrary constants, are general solutions of Eqs. and (155|) . respec- 
tively. At H//3 > 0, these functions tend to zero, for t — oo. Substituting these functions 
into the other equations, we get 

= uJ^e-*"'"^ + dae-"""', (60) 

ft.(i)^d.- «^°"-'^-f^-f''°^ -'«.-/.. ,61) 

6pHo{l + ^o) 

where ^2, d^, and (^4 are arbitrary constants. The two last expressions are valid for /3 7^ 4/3. 
For 13 = 4/3, 

\ 1 + ^0 / 

We see that none of the perturbations tends to infinity at t — )■ 00, so that the de Sitter 
solutions are stable. We should note that the fixed point, which corresponds to a de Sitter 
solution with fixed Hq, is not isolated, because there is only one condition: A = 3-?/q(1 + ^o) 
on the two arbitrary parameters Hq and ^o- Changing ^0, we get a new value of Hq for the 
same A. This is the reason why the function hi{t) has an arbitrary parameter d2- To analyze 
the stability of such a fixed point one cannot use Lyapunov's theorem |24j. 



10 



4.1.2 The Bianchi I metric 



The Bianchi universe models are spatially homogeneous anisotropic cosmological models. 
There are strong limits on anisotropic models from observations [25]. Anisotropic spa- 
tially homogeneous fluctuations have to be strongly suppressed, and models developing large 
anisotropies should be discarded as early or late times cosmological models. 

Interpreting the solutions of the Friedmann equations as isotropic solutions in the Bianchi 
I metric, we include anisotropic perturbations in our consideration. A similar stability 
analysis has been made for cosmological models with scalar fields and phantom scalar fields 
in [26]. The stability analysis is essentially simplified by a suitable choice of variables. Let 
us consider the Bianchi I metric 

ds"^ = - dt"^ + al{t)dxl + al{t)dxl + al{t)dxl. (63) 

It is convenient to express in terms of new variables a and Pi (we use the notation of [27]): 

ai{t) = a{t)e^^^'\ (64) 

Imposing the constraint 

/3i(t) + /32(t) + /33(t) = 0, (65) 
at any t, one has the following relations 

a(t) = (al(t)a2(t)a3(t))^/^ (66) 

H, = ^ = H + P,, and H = - = 1{H, + + H^). (67) 

Note that (3i are not components of a vector and, therefore, are not subject to the Einstein 
summation rule. In the case of the FLRW metric: ai = 02 = as = a, all Pi are equal to 
zero and H is the Hubble parameter; thus, we can use the same notations as for the FLRW 
metric. Following [27], we introduce the shear 

a' ^ Pi + Pi + Pi (68) 

It is easy to calculate that 

R = 12H^ + 6H + a'^, (69) 

Roo = -3 + if) - a^, (70) 

Rn = 9niHj + SHH,) = g,, (if + 3H' + (3, + 3HP,>j . (71) 
In the Bianchi I metric, we get □// = i? as 

ip= -?,Hi)-l2H'^ -QH -a'^. (72) 
Equations ([6]) and (fTOl) are now 

d = - md + ^ {l2H^ + 6H + (T^) , (73) 
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Pni = - -^Pm- (74) 

To get Eq. flTijl we have used condition (135|) on Wja- The Einstein equations have the form: 



a 



1 • 



:i + - + 2^v^ - 3^ - ^) + ^ + «^Vm = ^ 



(75) 



2H + 3H' + (3j - 3H(3j 



1 



2 '^^ 

Adding Eqs. ( !76|) for j = 1, 2, 3 and using (165|) . we get 



1 + - O + 2^V^+0-e+(2i/-/3,)(0-e) = A-/t'Pm. (76) 



a 



2H + 3H^ + — 



and adding now (!75i) and (1771) . 

2ij + a^] (1 + - + + - e - ^(0 - = -/^'(l + W^)pra. 

Assuming that P6|) is satisfied, we write Eq. (!78|) as 

1 



(77) 



(7J 



2 

24 



1 + 



+ 



4^2 



1.^0-.- 



(79) 



0" 



Subtracting (1761) . with j = i, from (176|) . with j = k, we obtain the following system of 
equations: 



A + 3H(3, -13k- 3H(3k 



:i + 0-O + (A-/3fc)(0-O = O 



^0) 



and using (!65|) . it is easy to get from this system 



(i + 0-O + A(^0^-^ 



and 



(l + 0-O + 2a^(-0^-^9 



0. 



f82) 



The functions H{t), o"^(t), 0(t), V'(^)) ^(^)) ''^(^) and Pm(^) can be obtained from equations 
fl39|l . (172D-([74]), ([79]) and ([82]). If H{t) and the scalar fields are known, then (3i can be directly 
obtained from ( ISTl) . 

The functions -ff(t), /3j(t), and cr'^it) are very suitable to analyze the stability of isotropic 
solutions in the Bianchi I metric. Indeed, the use of these variables makes the analysis of 
stability in the FLRW and Bianchi I metrics similar, because the equations of motion in the 
Bianchi I metric with = are identical to the equations of motion in the FLRW metric. 
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Thus, we can use some results from the previous subsection. In the neighborhood of the 
fixed point, which corresponds to de Sitter solution, we have the expansions ( 148|) - (!53|) . and 



(83) 



To first order in e, we get the following system of equations. Equations ( !54|) - (!56!) are valid 
for the Bianchi I metric as well, and instead of Eqs. ( !57l) and ( |58l) we get the system: 



d 
dt 



hi 

^1 



/3 



/3 



- 1 - - WJo>l - TT^P^ 



- SHotpi - 12Hohi 

- 6Hoal. 



0^2 



3/3 
12 



'ml 



cr 



;i+eo 



1) 



1 - 



(3 



K 



From the last equation, we get 



a 



-eHot 



(84) 
(85) 
(86) 

(87) 



where d^ is an arbitrary constant. 

The expressions for the functions (f)i{t), pi{t) and '&i{t) in the Bianchi I metric coincide 
with the corresponding expressions in the FLRW metric, which are given by fl5^ - fl5Ul) . 
Substituting these functions into Eqs. (l84l) and (!85l) . we obtain 



hi 



6/3iJo(l + eo 



5 ^-6Hot 



12H. 



^1 



2(/3 - 2){<6Hll3di - UH^di - K^f3do 



ffo/3(3/3-4)(l + eo) 
and, at /3 = 4/3, we get 

'{SH^di + K^do)t 



3i/n 



^9) 



^1 



+ d. 



5 g-6Hot_ 



The function ^i(t), which can be defined using Eq. (J75l) . is a bounded function if all other 
first-order corrections are bounded. Thus, we come to the conclusion that de Sitter solutions 
are stable if ifo > and /3 > 0. So, the stability conditions in the cases of the FLRW and 
Bianchi I metrics coincide. 



4.2 The case A = 0. Normalized variables 



To analyze the stability of the de Sitter solutions at A = 0, following [T^ we transform the 
system of equations using new dependent variables 

""—jh' ^'-m- "--'-IF' ^-#7 
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and the independent variable N: 



d _ d 1 d 

dN^"'d^^HJf ^ ' 



The use of the Hubble-normahzed variables [28j and as time variable makes the equation 
of motion dimensionless. Note that a change of dependent and independent variables of 
this kind is actively used in cosmological models with scalar fields, in order to analyze the 
stability in the FLRW metric [291 EO] , as well as in models of inflation (see pT] and references 
therein). Clearly, 

dX _ 1 dX _ I f i] H \ _ 7] XdH 

dN ~ H~dr ~ ~ m [h ~ J ~ ~ 4IP~ HdN' ^ ' 

^ = L(J iL\=^ + ixw-^^ m) 

dN H yQHf 6H^f QHpJ 6fH^ (3 H dN' ^ ' 

Equations (IT^ and (fT5]l are equivalent to the following ones, in terms of the new variables, 
dX , 1 /3 , \ dH , , 

'^"^ 2.1 + 2»'.Y)-3iy + ifi-»'ll (95) 



dN (3' ' H \(3 J dN 

and Eq. (ITUl) can be written as 

dZ A, , ZdH 

= -X-lZ-2 . 96 

dN ' HdN ^ ' 

To get the full system of equations we need one for In terms of the new variables, Eq. ([8]) 
has the form 

m ~ZY - AXW +^X + Z = 1. (97) 



Differentiating it 



usmg 



dW dY ^^^dX ^dW AdX dZ ^ 
'lN-'dN-''^dN-'''m+-pdN + dN='^ 

dY .f2XY_^ 



dN \ /3 
and substituting (|94l)-(l971). we get 



(4 fi \ 1 d FI 2 8 

AWX -6W--X + --Zj-— + 12W{X - 1) + -(6 - 4X - Z) - -X' = 0. (99) 

For Co = (and (3 ^ 4/3) the de Sitter solution in terms of new variables corresponds to the 
following fixed point: 

H = H„ iyo = -J_, Xo = l, Yo = ^^, Zo = ^^^j^. (100) 
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For P = 2, when po = 0, the stabihty of the de Sitter solution has been proven in [15]. In 
this paper, we discuss stabihty in the case /3 > 0. 

Let us consider perturbations in the neighborhood of the de Sitter solution (11000 : 

X = Xo(l + exi{N)) + 0{e^), Z = Z^{1 + ezi{N)) + 0{e^), 

W = Wo{l + ewi{N)) + 0{e^), H = Ho{l + sh{N)) + 0(6^). 
To first order in e, we obtain the system of linear equations: 

^=-3x, + l^, (101) 
dzi 4 dhi 

4^,,_l^^^+f|_3V... (103) 



dN 13 2(3 dN \f3 

djH^_ 8(/3-4) _ 2(3/3-4)(/3-2) 

dN /3(3/32-ll/3 + 12) ' /3(3/32-ll/3 + l2) ^ ^ 

Substituting ffTOlj) into ffTOTj) and (fT02D . we get 

dxi _ (/3-l)(3/3-4)2 _ (3/3-4)(/3-2) 

diV ~ " /3(3/32-ll/3 + 12)^' " /3(3/32_ii;3 + 12)^1' 

d^i _ 4(3/3^ - 7/3 - 4) 4(3/3-4)(/3-2) ^ ^ 

diV ~ /3(3/32-ll/3 + 12)'^' ^ /3(3/32-ll/3 + 12)^'' 

which has the following solution: 

a;, = cie^i^ + cse^^^, Zi = Cge^^^ + 046^^^, (106) 
where Ci and C2 are arbitrary constants 

_ _ 9/33 - 21/32 + IQ-D _ _ 9/33 - 21/32 + 16 + D 

~ 2(3/3-4)(/3-2) ~ 2(3/3-4)(/3-2) 

_ 9/33 - 45/32 + 80/3 - 48 + D \ _ _ 9/3^ - 45/32 + 80/3 - 48 - D 
^ ~ 2/3(3/32-11/3 + 12) ' ^ ~ 2/3(3/32-11/3 + 12) ' 



D = V8I/36 _ 378^5 + 297/34 + 1104/33 - 1984/32 + 256/3 + 768. 
We see that Ai = at 

/3i = 4 /32 = 2, /33,4 = ^±^z. (107) 
3 6 6 

The real part of Ai is negative in the interval /3 G (4/3, 2). It is easy to show that also the 
real part of A2 is negative in this interval. 

Therefore, the perturbations Xi and zi decrease in 4/3 < /3 < 2. Substituting a;i(A^) and 
zi{N) into Eqs. (11031) and (11041) . we get that hi{N) and Wi{N) decrease as well. Note that 
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hi{N) has a part Hi, which does not depend on N, and therefore can be considered as part 
of Hq. This result corresponds to the fact that, for A = 0, the value of Hq can be selected 
arbitrarily; thus, one can choose Ho = Ho + Hi instead of Hq. Adding here the results of 
|15] . we can summarize that the de Sitter solutions are stable with respect to perturbations 
in the FLRW metric with 4/3</3^2. At/o>0 the stable de Sitter solution corresponds 
to po ^ 0. 

Consider now the case of an arbitrary cq. For the de Sitter solution, we get 

H-H X -^ 7 - ^(^-^) 
n — no, Aq — i, ^0 — ; 

Y = -A ^g-{3~4//3)(7V-7Vo) u/ = _Jl ^g-{3~4//3)(7V-7Vo) 

3/3-4 9i/o/o ' 3/3-4 QHofo 

where No = Hoto- For /3 > 4/3, 

(3 2 

lim Y = — 7, lim W 



N^oo 3/3-4 N^oo 3/3-4 

Thus, the de Sitter solution tends to a fixed point, which means that, for any e > 0, there 
exists a such number Ni that the de Sitter solution is in the e/2 neighborhood of the fixed 
point for all N > Ni. Therefore, the stability of the fixed point guarantees the stability of 
de Sitter solutions for any value of Cq. We reach the conclusion that, at A = 0, all de Sitter 
solutions are stable for 4/3 < /3 ^ 2. 

For /3 = 4/3, we have no fixed point, because the Y and W corresponding to de Sitter 
solutions are proportional to N. Thus, this choice of dimensionless variable is not suitable 
to consider stability of the de Sitter solutions for /3 = 4/3. For /3 < 4/3 and /3 > 2, we find 
that the solutions are in fact unstable. 



5 Conclusions 

In this paper, we have investigated de Sitter solutions in a nonlocal gravity model, which 
is described by the action given in (see also [I5]). To carry out this task we have used 
the local formulation of the model ([3]), which includes two scalar fields. De Sitter solutions 
play a very important role in cosmological models, because both inflation and the late-time 
Universe acceleration can be described as a de Sitter solution with perturbations. 

We have found the interesting result that the model has de Sitter solutions only if /(r^) 
satisfies the second-order linear differential equation (fT6l) . If we consider models without 
matter or with a perfect matter fluid with a constant EoS parameter ^ ~l/3, then 
f{r)) can be an exponential function or a sum of exponential functions. For the model with 
f{rj) equal to a sum of exponential functions, particular de Sitter solutions have been found 

in ^5im\. 

In this paper, we have considered de Sitter solutions for an exponential /(r/) and found 
all solutions in the FLRW metric that correspond to a constant, nonzero value of the Hubble 
parameter, Hq. In particular, we have obtained the de Sitter solutions in the case when the 
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matter included in the model is dark matter. This case has never been considered before in 
the literature. 

When t — )■ oo, the de Sitter solutions tend to fixed points. In the model considered, the 
value of the cosmological constant does not fix Hq, therefore, the fixed points that correspond 
to de Sitter solutions are not isolated. We have analyzed the stability of these solutions in 
the cases of FLRW and Bianchi I metrics and obtained that the first-order corrections have 
no increasing mode^, this being valid for any nonzero value of A, and for Hq > and (3 > 0. 
They display constant and decreasing modes. For this reason, we can say that, for Hq > 
and /3 > 0, our solutions are stable for all nonzero values of A. For A = 0, the stable solutions 
correspond to i^o > and 4/3 < /3 ^ 2. 

Looking further, it will be interesting to consider the stability of the de Sitter solutions 
and the corresponding ghost-free conditions in the Einstein frame, for models with more 
general functions /(r^) satisfying the differential equation (fT6|l . 
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